Precise Measurement of the 7r+ 7r°e+i/ Branching Ratio 
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Using a large acceptance calorimeter and a stopped pion beam we have made a precise mea- 
surement of the rare tt^ — > Tv'^e'^v (tt^) decay branching ratio. We have evaluated the branching 
ratio by normalizing the number of observed irp decays to the number of observed tt^ — > e^i^ (''"e2) 
decays. We find the value of r(7r+ 7r''e+i/)/r(total) = [1.036 ± 0.004 (stat.) ± 0.004 (syst.) ± 
0.003 (7re2)] X 10~*, where the first uncertainty is statistical, the second systematic, and the third is 
the 7re2 branching ratio uncertainty. Our result agrees well with the Standard Model prediction. 

PACS numbers: 13.20.Cz, 11.40.-q,14.40.Aq 



The rare pion beta decay, tt"^ — > T:^e^v (branching ra- 
tio R-^p ~ 1 X 10~^), is one of the most basic semileptonic 
electroweak processes. It is a pure vector transition be- 
tween two spin-zero members of an isospin triplet, and 
is therefore analogous to superallowed Fermi (SF) tran- 
sitions in nuclear beta decay. Due to its simpHcity and 
accuracy, the theory of Fermi beta decays is one of the 
most precise components of the Standard Model (SM) of 
electroweak interactions. 

The CVC hypothesis P, Q and quark-lepton universa- 
lity relate the rate of pure vector beta decay (both pion 
and nuclear) to that of muon decay via the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix ele- 
ment Vud Including loop corrections, 5, the rate of 
pion beta decay is given by 0, 0| : 
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where is the Fermi weak coupling constant, A — 
Af+ — A/o, e — (me/A)^, Af+, Mq, and mg are the masses 
of the 7r+, 7r°, and the electron, respectively, while /, the 
Fermi function, is given by 
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The main experimental source of uncertainty in P^^ 
amounts to just 0.05%; it comes from the measurement 
of A 0] ■ The combined radiative and short-range physics 
corrections amount to (5 ~ 0.033 and are exceptionally 



well controlled, yielding an overall theoretical uncertainty 

of of < 0.1% UlElili^. Hence, pion beta decay 

presents an excellent means for a precise experimental 
determination of the CKM matrix element Vudi hindered 
only by the low branching ratio of the decay. 

The CKM quark mixing matrix has a special signifi- 
cance in modern physics as a cornerstone of a unified 
description of the weak interactions of mesons, baryons 
and nuclei. In a universe with three quark generations 
the 3x3 CKM matrix must be unitary, barring certain 
classes of hitherto undiscovered processes not contained 
in the Standard Model. Thus, an accurate experimental 
evaluation of the CKM matrix unitarity provides an in- 
dependent check of possible deviations from the SM. As 
the best studied element of the CKM matrix, Vud plays 
an important role in all tests of its unitarity. However, 
evaluations of Vud from neutron decay have, for the most 
part, not been consistent with results from nuclear SF 
decays 1 1 . Clearly, a precise evaluation of Vud from 
pion beta decay, the theoretically cleanest choice, is of 
interest. 

The most precise measurement of the pion beta de- 
cay rate on record was made by McFarlane et al., at 
LAMPF by detecting in-flight 7r+ decays in the LAMPF 
7r° spectrometer This work reported Ptt/s = 0.394 ± 
0.015 s""'^, which is an order of magnitude less precise 
than the theoretical description of the same process. 
Hence, we initiated the PIBETA experiment, a program 
of precise measurements of the rare pion and muon de- 
cays at rest, chief among them being the pion beta decay, 
at the Paul Scherrer Institute (PSI), Switzerland 0|. 

In this Letter we present an analysis of the tt"^ — > 
n^e^v decay events recorded with the PIBETA appara- 
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FIG. 1: A schematic cross section of the PIBETA detector 
system. Symbols denote: BC-thin upstream beam counter, 
ACl,2-active beam collimators, AD-active degrader, AT- 
active target, MWPCl,2-thin cylindrical wire chambers, PV- 
thin 20-segment plastic scintillator barrel. BC, ACl, AC2, 
AD and AT detectors are also made of plastic scintillator. 



minimum-bias sets, we extract the tt^ and event sets, 
the latter for branching ratio normalization. In a stopped 
pion experiment these two channels have nearly the same 
detector acceptance, and have much of the systematics in 
common. 

A full complement of twelve fast analog triggers com- 
prising all relevant logic combinations of one- or two-arm, 
low- or high calorimeter threshold (labeled ht and lt, re- 
spectively), prompt and delayed (with respect to Tr~^ stop 
time), as well as a random and a three-arm trigger, were 
implemented in order to obtain maximally comprehen- 
sive and unbiased data samples. 

Signal definition and accurate counting of the TTe2 
events for normalization present a major challenge in this 
work. As in all previous studies, our TTe2 data include 
undiscriminated soft-photon 7Te2-y events. Due to positron 
energy straggling in the target, accidental coincidences of 
multiple muon decay events, and the calorimeter energy 
resolution function, the Tre2 events are superimposed on a 
non-negligible muon decay background. This background 
was removed by fitting the measured e+ timing spec- 
tra with the functions for pion decay (signal), muon de- 
cay (background) , plus the associated pile-ups (see Fig. |3 
top) . We also extracted the absolute Tre2 branching ratio 
using this method and normalizing to the number of pion 
stops in the target. The results were in agreement with 



tus from 1999 to 2001. We tuned the ttEI beam hue 
at PSI to deliver - 10^ 7r+/s with ~ 113MeV/c. 
The pions were slowed in an active degrader detector 
(AD) and stopped in a segmented 9-element active tar- 
get (AT), both made of plastic scintillator. The major 
detector systems are shown in a schematic drawing in 
Fig. n Energetic charged decay products are tracked in 
a pair of thin concentric multiwire proportional cham- 
bers (MWPCs) and a thin 20-scgment plastic scintillator 
barrel veto detector (P V) . Both neutral and charged par- 
ticles deposit most (or all) of their energy in a spherical 
electromagnetic shower calorimeter consisting of 240 ele- 
ments made of pure Csl. The entire detector system, its 
response to positrons, photons and protons, energy and 
time resolution, signal definitions, along with other rel- 
evant details of our experimental method, are described 
at length in Ref. 

The measurement relies on detecting the 7r° 77 de- 
cay which immediately follows a pion beta decay event. 
The two photons are emitted nearly back-to-back, with 
about 67MeV each. Thus, the experiment is set to 
record all large-energy (above the /i — *■ evD endpoint) 
electromagnetic shower pairs occurring in opposing de- 
tector hemispheres during a ~ 180 ns long "pion gate" , 
ttG (non-prompt two-arm events). The ttG is timed so as 
to include a sample of pile-up events preceding the pion 
stop. In addition, we record a large prescaled sample of 
non-prompt single shower (one-arm) events. Using these 
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FIG. 2: Top panel: A typical histogram of time differences 
between the beam pion stop, f(DEG), and 1-arm ht event 
time, t{e), (dots), compared with a sum of the Monte Carlo- 
simulated responses for decay (tt), muon decay (fj.), and 
muon pile-up events {fJ.p). The ■Ke2 pile-up background, being 
much lower, is off scale in the plot. Prompt events are sup- 
pressed. Bottom panel: Csl calorimeter energy spectrum for 
the 7re2 decay events, after background subtraction. 
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FIG. 3: Histogram of time differences between the beam pion 
stop and the tt/j decay events (dots); curve: pion hfetime. 



FIG. 5: Histogram of the 7-7 opening angle in tt^ decay. 



the recommended Particle Data Group (PDG) value [ij 
at a sub-percent level, with the uncertainty dominated 
by the systematics of the stopped pion counting. The 
latter is absent in our determination of R-„f3. The 'Ke2 
energy spectrum after background subtraction is given 
in Fig. 13 bottom. The statistical uncertainty of the ex- 
tracted number of 1:^2 events, A^7re2, is negligible. 

The TT^ signal definition was more straightforward, as 
seen in Figs. 01 and ^ which show the pion decay time 
spectrum and 7-7 relative timing histogram, respectively, 
for TT^ events, both free of backgrounds. Finally, the 
histogram of recorded 7-7 opening angles for pion beta 
events, shown in Fig. [SJ provides a sensitive test of the 
systematics related to the geometry of the beam pion 
stopping distribution, an important contributor to the 
acceptance uncertainty. 

The Tip branching ratio R-^fs was evaluated from 
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where N-^p is the number of detected tt^ events, Nj^+ is 
the number of the decaying tt+'s, /ttg is the delayed pion 
gate fraction, A^J is the ht detector acceptance eval- 
uated by GEANT simulation, tj is the detector live time, 
/cpp is the correction due to the charged particle (CP) 
veto system pile-up, /d is the tt'^ Dalitz decay correction, 
and /ph is the photonuclear absorption correction. 
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FIG. 4: Histogram of 7-7 time differences for tt/j decay events 
(dots); curve: fit with a Gaussian function plus a constant. 



ev branching ratio i?Tre2 is given by 



R-Ke2 — 
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where p7re2 is the prescaling factor applied to 7re2 trig- 
gers, A^]2 is the high-threshold detector acceptance for 
TT — > ej^ decay events, including radiative corrections, 
while epv, ecij and ec2 denote the plastic veto and wire 
chamber efficiencies, respectively. Clearly, taking the ra- 
tio i?7r/3/^7rei/ Icads to Cancellations of many common 
factors, apart from small corrections taking into account 
slight differences in thresholds, trigger timing (two-arm 
vs. one-arm), weighting of the efficiencies, and similar ef- 
fects. Most importantly, N^+, the number of stopped 
pions drops out. The main sources of uncertainty are 
listed with their values in Table HI 

As the external systematic uncertainties are self-ex- 
planatory, we turn to the internal ones. The syste- 
matic uncertainty in NT^e2 comes from the muon-decay 
background subtraction discussed above, and reflects the 
propagated error limits of the method. The precision of 
A^J M^e2 is dominated by the uncertainty of the x-y-z 
distribution of pion stops in the target. The latter was 
determined with better than 50 /xm accuracy by tomo- 



TABLE I: Summary of the main sources of uncertainty in the 
extraction of R^p, given in % (see text for discussion). 



Type Quantity 


Value 


Unc. (%) 


external: r(7re2) 

r(^» -> 77) 

tt"'" lifetime 


1.230 X 10"" 
0.9880 
26.033 s 


0.33 
0.03 
0.02 


combined external uncertainties: 




0.33 


internal: A.^e2 systematic 

/IHT /4HT 

■Kp accid. bgd. 
/cpp correction 
/ph correction 


6.779 X 10** 
0.9432 
1.130 
0.00 
0.9951 
0.9980 


0.19 
0.12 
0.26 
< 0.1 
0.10 
0.10 


combined internal uncertainties: 




0.38 


statistical: N-ni3 


64,047 


0.395 
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graphic back-tracing of 7re2 and muon decay positrons 
into the target 0. Corrections due to the undetected 
low portions of the e and 7 energy spectra in the calorime- 
ter (the energy "tail") contribute weakly to the accep- 
tance uncertainty due to strong correlations between the 
energy responses to the two decay channels. This ex- 
periment has a unique advantage over its predecessors: 
it measures branching ratios as well as differential an- 
gular and energy distributions of decay products for all 
rare pion and muon decays simultaneously. This provides 
multiple redundant consistency checks of the evaluated 
and simulated acceptances (cf., e.g., Ref. (3)- the 
present analysis the largest internal contribution to the 
systematic uncertainty comes from the ratio of gate frac- 
tions, r^rG , due to our decision to include even the earliest 
■Kp decay events, thus maximizing the useful event statis- 
tics. An alternative method involves a sharp cut at, say, 
10 ns after the pion stop time, which reduces the syste- 
matic error, along with the tt^ event statistics, yielding, 
however, a consistent R-rrfi- The inherent resolution in 
the zero time point is excellent — it relies on the prompt 
A(7r+, 7r*')B signal and the accelerator rf pulse, providing 
timing calibration at the level of ~20ps or better, and 
room for further improvement of the r^^Q precision. The 
pile-up correction /cpp was evaluated using a random 
trigger, and confirmed by simulations. We modified our 
GEANT3 code to calculate the photonuclear correction 
/ph, and conservatively assigned it a 50% uncertainty 
(details are given in Ref. pjj). Efficiencies epv, ecij and 
ec2, not listed in Tabled were measured with an accu- 
racy of 0.01% fiij. 

Using the above method and the PDG 2004 recom- 



mended value of i?° e2 



1.230(4) X 10- 



we extract 



the pion beta decay branching ratio, our main result: 

R°^^ = [1.036 ± 0.004 (stat) ± 0.005 (syst)] x 10"* , (5) 
or, in terms of the decay rate. 



= [0.3980 ± 0.0015(stat) ± 0.0019(syst)] s" 



(6) 



which represents a six-fold improvement in accuracy over 
the previous measurement ,,12J . Alternatively, the nor- 
malization can be tied to the theoretical value i?^g2 = 
(1.2352 ±0.0005) x 10"'' ^ which would increase the ex- 
tracted by 0.42 % to 1.040 x 10"^ In a direct evalu- 
ation of the pion beta decay branching ratio using Eq ij^J, 
i.e., normalizing to the number of beam pion stops, we 
obtain R^p ■ 10* = 1.042 ± 0.004(stat) ± O.OlO(syst), con- 
sistent with our main result given in Eq. 

Whether scaled to the experimental or theoretical 
i?^e2, our result for E^^^ is in excellent agreement with 
predictions of the SM and CVC given the PDG recom- 
mended value range for Vud (ill : 



and represents the most accurate test of CVC and 
Cabibbo universality in a meson to date. Our result 
confirms the validity of the radiative corrections for the 
process at the level of 4(Texp, since, excluding loop cor- 
rections, the SM would predict Riy""^- = (1.005 - 
1.007) X 10-* at 90% C.L. 



Using our result, Eq. jSJ, we can calculate a new value 
of Vud from pion beta decay, l^^J^^^™^ = 0.9728(30), 
which is in excellent agreement with the PDG 2004 aver- 
age, V^J™*^ "''^ = 0.9738(5). We will continue to improve 
the overall accuracy of the tt^ decay branching ratio to 
^ 0.5 % by further refining the experiment simulation and 
analysis, and by adding new data. 
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